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Abstract —Cell association scheme determines which hase 
station (BS) and mobile user (MU) should he associated with 
and also plays a significant role in determining the average 
data rate a MU can achieve in heterogeneous networks. 
However, the explosion of digital devices and the scarcity 
of spectra collectively force us to carefully re-design cell 
association scheme which was kind of taken for granted 
before. To address this, we develop a new cell association 
scheme in heterogeneous networks based on joint consideration 
of the signal-to-lnterference-plus-nolse ratio (SINK) which a 
MU experiences and the traffic load of candidate BSfl. MUs 
and BSs in each tier are modeled as several independent 
Poisson point processes (PPPs) and all channels experience 
independently and Identically distributed {i.i.d.) Rayleigh fad¬ 
ing. Data rate ratio and traffic load ratio distributions are 
derived to obtain the tier association probability and the 
average ergodic MU data rate. Through numerical results. We 
find that our proposed cell association scheme outperforms 
cell range expansion (CRE) association scheme. Moreover, 
results indicate that allocating small sized and high-density 
BSs will Improve spectral efficiency if using our proposed cell 
association scheme in heterogeneous networks. 

Index Terms —Heterogeneous networks; cell association 
scheme; traffic load; Poisson point processes 

I. Introduction 

Driven by a new revolution of digital devices like smart 
phones, tablets and so on, there has been experiencing 
a tremendous growth of mobile internet traffic in recent 
years. Traditional network expansion techniques like cell 
splitting are often utilized by telecom operators to achieve 
the expected throughput, which are less efficient and proven 
not to keep up with the pace of traffic proliferation in 
the near future. Heterogeneous networks then become a 
promising and attractive network architecture to settle this. 
Heterogeneous networks are a broad term that refers to the 
coexistence of different networks (e.g., traditional macrocell 
and small-cell networks like femtocells and picocells), each 
of them constituting a network tier. Due to differences 
in deployment, base stations (BSs) in different tiers may 
have or use different transmission power levels, radio ac¬ 
cess technologies, fading environments and spatial densities. 
Heterogeneous networks are envisioned to cope with most 
problems of existing network architectures like dead spots, 
inter-cell interference, less efficient, etc and has been intro¬ 
duced in the LTE-Advanced standardization [I]. Massive 

'Candidate BSs is comprehended as the set of BSs which a mobile user 
(MU) is most likely associated with. 


work has been done in heterogeneous networks scenario 
mainly related with coverage modeling El, El, cooperative 
communications H, energy consumption modeling El, 0, 
interference cancellation Q, interference management ||8l 
and resource allocation a-iini, however none of which 
pays enough attention to existing problems on cell associa¬ 
tion schemes. 

A. Motivation and related work 

In heterogeneous cellular networks, there are more BSs 
which a MU can choose to be associated with than in tradi¬ 
tional homogeneous single-tier cellular networks. Therefore, 
cell association scheme is an indispensable factor in wireless 
networks modeling. By using the maximum received signal 
strength (RSS) as cell association scheme, ElSawy and 
Hossain quantified the performance gain in the outage prob¬ 
ability obtained by introducing cognition into femtocells in 
two-tier heterogeneous networks ||2l. In IfT^ . Ali and Saquib 
developed a practical yet tractable method of evaluating 
vertical handover algorithms in a WLAN/Cellular two-tier 
heterogeneous network and the cell association is also based 
on the maximum RSS. Dhillon et al. El proposed a tractable 
and accurate model for a downlink heterogeneous cellular 
network consisting of K tiers of randomly located BSs. 
Novlan et al. aimed to evaluate two fractional frequency 
reuse (EER) methods - strict EER and soft frequency reuse 
by using Poisson point processes (PPPs) in llT3l . The cell 
association scheme utilized by Dhillon and Novlan is based 
on the maximum downlink signal-to-interference-plus-noise 
ratio (SINR). Also, the nearest BS cell association scheme 
is applied in some literatures like El, mi. Yong Sheng 
et al. investigated the design and the associated tradeoffs of 
energy efficient heterogeneous cellular networks through the 
deployment of sleeping strategies in ||5l- In IfT^ . Mukherjee 
provided a general theoretical analysis of the distribution 
of the SINR at an arbitrarily-located MU in heterogeneous 
networks. 

Prom above literatures, existing cell association schemes 
have been mainly based on the RSS, SINR or the distance 
from nearby BSs to determine which BS and MU should be 
connected with each other. This is legitimate for traditional 
homogeneous single-tier cellular networks where the RSS 
or the SINR serves as a good indicator of the data rate 
received by the MU. However, it is no longer the case 
in heterogeneous networks in which BSs from different 
tiers transmit wireless signals at very different power levels, 


varying from milliWatt (mW) to Watt (W); a) the higher 
RSS may be a result of the higher transmission power 
used by the BS. It may cause congestions in BSs which 
have higher RSS and idleness in BSs whose RSS is lower 
whereas can still guarantee successful transmission. This 
result brings unbalance and inequity among BSs in different 
tier networks; b) the number of MUs served by a small¬ 
cell BS is typically small due to its much smaller coverage. 
Consequently, the current traffic load of the BS plays a 
significant role in determining the share of BS capacity 
received by each MU. For example, the joining of a MU 
into a small-cell BS currently serving one MU may halve 
the data rate received by the current MU; c) as for choosing 
the nearest BS for association, it is so impractical that only 
used for theoretical analysis. Thus, it is no longer optimum 
to determine cell association solely based on the RSS, the 
SINR or the distance from nearby BSs. 

As described above, cell association schemes play an 
important role in determining the allocation of spectral 
resource in BSs, the transmission rate that a MU can achieve 
and even the energy consumption of MUs. In ca, m, 
a solution was proposed to partially solve the problem 
a) by introducing a biased factor or B into the RSS, 
which allows an expansion of the coverage of small-cell 
BSs. The effectiveness of the scheme however remains 
questionable in networks with inhomogeneous user density, 
e.g. MUs clustering around BSs. In iflTll . authors mentioned 
the problem b) in the subsection of resource allocation. 
However, per MU data rate is only a performance metric 
with the form of rate coverage and the used cell association 
scheme was still conventional, which left these problems 
unsolved. 

B. Contributions and organization 

To solve problems a), b) and c), a spectrum efficient 
cell association scheme based on the joint consideration of 
the received SINR and the traffic load of BSs is proposed 
for heterogeneous networks. To match real BSs deployment 
scenarios, PPP is used to model heterogeneous cellular 
networks in this article, which has been strengthened by the 
empirical validation IfTSl and the theoretical validation lfT9ll . 
The contributions and novelties of this paper are summarized 
as follows. 

1) A new cell association scheme is proposed with two 
steps for heterogeneous networks. The first step is 
mainly for choosing the candidate BSs by traditional 
method, i.e., the nearest n BSs, while the second 
step determines the ultimate one BS based on the 
consideration of the received SINR experienced by a 
MU and the traffic load of candidate BSs. 

2) Following the cell association scheme and taking a 
three-tier heterogeneous network as an example, the 
tier association probability and the average ergodic 
MU data rate are derived for numerical analysis. 

3) Based on numerical results, the new cell association 
scheme outperforms CRE association scheme. 

The reminder of this paper is organized as follows. In 
section m we present our network model and propose a 


new cell association scheme in general case. A three-tier 
heterogeneous network is analyzed in section [HIl Moreover, 
the tier association probability of heterogeneous networks 
is derived for performance analysis. Section |IV] presents the 
numerical results of the proposed cell association scheme. 
Section [Y1 concludes this paper. 

II. Network model and proposed cell 

ASSOCIATION SCHEME 

We Consider a K-tier heterogeneous downlink cellular 
network which consists of macrocells, picocells, femtocells, 
etc. BSs of each tier are assumed to be spatially distributed 
following independent homogenous PPPs denoted by <1)^, 
k G - ,iT}. The BS intensity of the fc-th tier network 

is Afe, k G {1, 2, • • ■ ,K}. MUs are located according to a 
homogeneous point process denoted by <i)ii with intensity 
A„. All BSs in the same tier network are configured with 
the same transmission power P^, k G {1,2, ,K} and 

share the same bandwidth. BSs in different tier networks 
are configured with different bandwidths. Moreover, within 
a cell, MSs are allocated by orthogonal frequencies. There¬ 
fore, there is no intra-cell interference in a cell. Also for 
simplicity, the open access policy is applied for MUs. It 
means all MUs can be served by BSs in any tier networks. 

We propose a new cell association scheme that bases its 
cell association decision on the instant traffic load of each 
BS and the transmission rate that can be allocated by the 
BS. More specifically, the cell association scheme can be 
divided into two steps. 

1) If a MU wants to be associated with a BS, it will 
firstly choose n nearest BSs from each tier as the 
candidate BSs. The candidate BS set is defined as 

= {(fc, i) |fc G (1, 2, • • • , AT} , i G [1, n]}, where 
k is the fc-th tier network in a heterogeneous network 
and i is the Uth BS in the n nearest BSs from the 
/c-th tier network. For example, (1,3) represents the 
3rd BS in the 1st tier network. The total number of 
BSs in Ub is nK. 

2) The MU will select a candidate BS from CIb- 
This selected BS will send data to the MU 
with the maximum average transmission rate, i.e., 

■ In (1 + SINRfe^i), where Bk,i and Nk,i are the 
total bandwidth of a candidate BS in set Hb and the 
instant BS traffic load, respectively. SINR^ ^ is the 
MU instant SINR associated with a BS {k,i) in the 
set Ub. It is assumed that the total BS bandwidth is 
shared equally among all associated MUs0 

III. Performance analysis 

Without generality, in the following we will analyse the 
scenario when K = 3 and n — 1. The network being 
considered is a three-tier heterogeneous network and only 
the nearest BS at each tier from a MU can be chosen 
as candidate BSs. The analysed three-tier heterogeneous 
network is depicted in Fig. [T] 

^In our following analysis, the candidate BS {k,i) can be denoted by 
k when n = 1. However, weTl keep using {k,i) for completeness and 
preciseness. 
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Figure 1. The three-tier heterogeneous network (20kmx20km ) 
modeled as a superposition of three independent Poisson Voronoi 
tessellations. These polygons are 1st tier cells (edges with red solid 
lines), 2nd tier cells (edges with blue dot-dash lines) and 3rd tier 
cells (edges with black dotted lines). 


A. The downlink SINK distribution 


/siNRfc.i (x) = 27rAfc r irr'^XkT' {x) + ^ 

exp [—Trr^Afc (1 + if (a;))] x 
exp (-^^) dr 


(3) 


with 



X 


1 

1 + a; 


(4) 


B. The tier association probability 

The spatial average data rate a MU can achieve is denoted 
by Ck, k G {1,2,3} and the event that the considered 
MU is associated with a BS in the fc-th tier network is 
denoted by {Nrier — k). We’ll apply Slivnyak’s theorem 
to the following analysis on a MU that located at the origin 
o, which implies that conditioning on having that user at 
the origin, properties of all coexisting PPPs maintaining the 
same. Using the above association scheme in section |II] the 
probability that a MU is associated with the 1st tier BS (1, i) 
is 


For downlink transmission of a BS (fc, i) to a MU which 
is located at the origin o, the SINK experienced by this MU 
is expressed by 


SINRfe,i 


Pkhk,ir 


— a. 

k,i 


E ^ Pkhk^mTk^ + fT2 ’ 


( 1 ) 


where is set of interferers in the fc-th tier network. 
hk,i and hk^m are channel power gains due to small-scale 
fading between the considered MU and BS {k,i), {k,m), 
respectively. For convenience and without generality, we 
assume hk,i ~ exp (1) and hfc.m ~ exp (1). The background 
noise is assumed to be additive white Gaussian with variance 
cr^. and are path losses with a being the path 
loss exponent, ^ and rk^m being the respective Euclidean 
distance to the corresponding BS (fc,z) and {k,m), respec¬ 
tively. 

Referring to 0, the coverage probability that a MU is 
covered by its nearest BS in a particular tier fc is derived as 
follows 


P^ (x) = Pr (SINRfc,, > :e) 

/■-too 

= 27rAfc / r exp [-Trr^Afc (1 + (a:))] X (2) 

Jr^O 


exp 




Pk 


dr, 


where p (x) = x^ — ^-nrdy. P^ (a;) is the complemen¬ 

tary cumulative distribution function (CCDF) of SINR^ ^ 
and the cumulative distribution function (CDF) of SINR^ ^ 
is 1 — P^ (a;). By taking a derivative of the CDF of SINR^^i 
with respect to x, the probability density function (PDF) of 
SINRfc i is obtained by 


Tl =Pr(7VT,er = l) 


= Pr 
= Pr[ 


Cl > maxCj 

In (1 + SINRi,i) > 


B2,i 
N2 , i + 1 


Wi.i-i-i 






ln(l-bSINR 2 ,.), 
ln(l-TSINR 3 ..)] 


ATl.i-l-l 


In (1 + SlNRi,i) 


Pr[i^ln(l + SlNRi.,)> 


a=2 

B. 


JVi.i + l 

In (1-b SINRj- 


Ni 


< ^ 
i + 1 ^ B, 


ln(l-|-SINRi,i) 

In(l-l-SINRj’i) 


N. 

3“ 

= n Pr 

= n Pr K, <|^-S 1 NRi/,- 

i=2 L 


n fo 

j=2 


+ CXD 




N- 


i/i 


(^■^) '-^SINRi/,- (x)dx 


(5) 

where (I) is due to the independence between two events 
> C^j and {Cl > ^}; in (II), Ni/^, SINRi/^ denote 
and respectively; (III) is obtained by 

applying the law of total probability where FaTj ,. (x) is 


the CDF of 




and /siNRw (x) denotes the PDF of 


Nj^i + l aoiiNKi/j 

ln(l-|-SINRi,i) 

In(l-l-SINRj’i) ■ 

1) The CDF of In this paper, it is assumed that 

each BS has a unique saturated downlink transmission queue 
for each MU. This assumption implies that MU always has 
data to receive from a BS which covers and associates that 
meH Each MU choose the associating BS with probabilities 
denoted by Tfc, fc € {1,2,3}. Thus, three point processes 
denoted by and are formed by thinning the 

original PPP The thinned processes are the locations 
of MUs which are associated with the 1st, the 2nd and the 


^The shape of a BS’s coverage is Voronoi-tessellated. 


































3rd tier BSs. The thinned point processes are still PPPs and 
intensities are 71 A„, 72 A„ and TsA^, respectively. 

Through interpretations above, the probability mass func¬ 
tion (PMF) of is given by 


fNk,i (n-) = Pr {Nk,i = n) 

^ r-too 

Jo n! 


fSk (s) ds 




r(c) 


-ds 


= fo 
_ mxur {cXkV 

“ n\ r(c) ^ 

j+oo „n+e-l^-s(rt,A„-rcAt.)^„ 

^ (rfcA„)"(cAfc)° r(n-rc) 

’^!r(c) (rfcA„-rcAfc)”+‘= ^ 

f+°° (TfcA„+cAfc)'‘+° —s(7fcA„+cAfc) j 

Jo r(n+c) 


2j The PDF of ' Let the PDF of SINR^,, be 

/siNRfc.i {.x) and the CDF be Fsinr^ ^ (x), then the CDF of 
In (1 -f SINRfe^i) is derived by 

^ln(l+SINRfc,i) {y) =Pr[ln(l-f SINRfe,,) < y] 

= Pr (SINRfe,i < et'- 1) . (10) 

= PsiNRfc,, (e^ — 1) 

By taking a derivative with respect to y in both sides of 
([Tol l, the PDF of In (1 -f SINR^ i) is obtained by 


/ln(l-|-SINRfc,i) {y) = e^- /siNRfc,i - 1) . (11) 


Let fjo (x, y) denote the joint probability density 
function (JPDF) of random variable tuple 
(ln(l + SINRi,i),ln(l + SINRj-i)), j G {2,3}. The 
in(i+siNR,..) L derived by 

/siNRi/,- {z) = \y\- fjo {zy, y) dy 


PDF of 


= /o°° y fjo {zy, y) dy 


( 12 ) 


0) (rfcA„)"(cAC‘ r(n+c) 

{Tk\,.+c\kF+‘ r(n+i)r(c) 

where P (•) is the Gamma function; fst, (s) 


(6) 


r(c 


is the PDF of the Voronoi cell area 


of the /c-th tier network obtained through simulations and 
c = 3.575 is a constant ll20l . (I) is obtained due to the 
integration of the formula with underline is 1 over the 
domain. Actually, the format of the formula with underline 

fc — 1 —^/9 

is the PDF of Gamma distribution like y = ^ • 


The CDF of 


Afi.i + l 

Ni,i + 1 


is derived by 


F, 


N, 


1/i 


{x) =Pr(^ 


,+i 


< X 


+ 1 

= Pr [Ni^i < (iVj-, + 1) cc - 1] 

m °° 

= E Pr [A^i,. < (f + 1) X - 1 I , = f] X 
t=o 

Pr(iV,-,=f) 

OO 

^Ni,i [L(^ + l)a; - IJ] •/7v,,i (t) 

t=o 


where ^ (1) is the CDF of Nk,i which is derived by 


(7) 


PNk,i (0 = fNk,i (n) 


n—0 


_ (TkAuFjcAk) L (n+c) j rp. N 

en (AA„+cAfc)’*+‘= r(n+l)r(c) > ^ t [U, OO) 


r(n-rc) 


n—0 


( 8 ) 

[• J is the floor function and (I) follows the law of total 
probability. Substituting (l8]l back in O, we obtain the CDF 
of as follows 


Ni/i 2^ l(7^.a„+cAj)‘+" r(t+i)r(c) ^ 

L(i+iV-ij ■ (9) 

vG (riA„)"(cAir r(K+c) 1 

(riA„-rcAi)’*+‘= r(n-ri)r(c) / 


r(t-re) 


^ in this paper, traffic load of a BS is defined as the total number of 
MUs associated with that BS. 


where (I) is obtained by using the ratio distribution (or 
quotient distribution) formula of two nonegative random 
variables. Because of the independence of the two variables, 
i.e.. In (1 -f SINRi^i) and In (1 + SINR^^i), fjo (x, y) is ex¬ 
pressed by 


fjo {x, y) — /ln(l-|-SlNRi i) (x) • /ln(l-|-SlNRj i) (l/) 

= e-+yfsmR ,. (e" - 1) • /sinr, . (e« - 1) ’ 

(13) 

Substituting (fTsT i back into (fT^ . we obtain the PDF of 

ln(l-|-SINRi,i) 

In(l-l-SlNRj’i) 


as follows 


/siNRi/j {z) = /p°° yel^+Lt//siNR - 1) x 

where fsiNRk,, {x) is given by 

Similarly, by repeating the derivations above, the proba¬ 
bilities that a MU is associated with a BS in the 2nd and 
the 3rd tier network are obtained, i.e.. 


Ts 


Ta 


= n /o ■fsiNRy,{x)dx , 

i=ij¥2 ^ ^ 

(15) 


=n /c 

i=i 


+ 00 


F 


'^3/. '/siNRa/, {y:)dx . (16) 


Using numerical method, the exact value of 7/c, k G {1,2,3} 
are obtained by solving equations ds), (fTSl l and (fThl l. 


C. The average ergodic MU data rate 

In this subsection, we focus on the average ergodic MU 
data rate of a MU in a K-tier heterogeneous network. We 
assume that Shannon’s capacity can be achieved by some 
coding methods. The average ergodic MU data rate can be 
obtained by considering per tier user data rates weighted by 
the corresponding tier association probabilities. The average 
ergodic MU rate in a 3-tier heterogeneous network is given 
by 

3 

K FA- 


(17) 
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Figure 2. The average ergodic MU data rate with respect to the 
density of the 2nd tier BSs with spectral allocation of (B1>B2>B3) 
in a three-tier heterogeneous network. 


Figure 3. The average ergodic MU data rate with respect to path 
loss exponent with two kinds of spectral allocations in a three-tier 
heterogeneous network. 


Kfc is the ergodic MU data rate conditioning on a MU is 
associated with a specific BS in the fc-th tier network which 
is given by 


Jifc =E[Bfc,,-ln(l + SINRfe,,)] 

= r Pr [In (1 + SINRfe,0 > t] dt , (18) 

= ^mCPc (e*-l)dt 

where (I) is derived because In (1 -|- SINR^ ^) is a nonneg¬ 
ative random variable; E (•) is an expectation operator and 
Pc (■) given by (|2]i. Substituting (fTSl l into (fTTl i. we can get 
the unconditional average ergodic MU data rate as follows 


3? 


exp [—Trr^Afe {1 + ip{e* - 


exp 




Pk 


drdt 


rx 

- 1))] X . 


(19) 


IV. Numerical results and discussions 


This section presents numerical results of previous sec¬ 
tions, followed by discussions. Parameters used in this 
article are refereed to existing work focused on hetero¬ 
geneous networks. Specifically, we assume that cr^ = 0 
which denotes a interference-limited scenario. BS densities 
and BS transmission powers are A 2 = 2Ai, A 3 = 20Ai, 
A„ = 50Ai, Pi = 53dBm, P 2 = 33dBm, P 3 = 23dBm US), 
Ha. Allocations of spectra are divided into 4 cases, i.e., 
(B1>B2>B3): Bi = 15MHz, B 2 = lOMHz, B 3 = 5MHz; 
(B1>B3>B2): Bi = 15MHz, B 2 = 5MHz, B 3 = lOMHz; 
(B2>B3>B1): Bi = 5MHz, B 2 = 15MHz, B 3 = lOMHz; 
and (B3>B2>B1): Bi = 5MHz, B 2 = lOMHz, B 3 = 
15MHz ITtII^ In the following, we will use default values 
above unless otherwise declared. 

Fig. la shows the average ergodic MU data rate with 
respect to the density of the 2nd tier BSs A 2 which varies 
from 0.4(fcm^)“^ to 0.9(fcTO^)“^ considering three different 
path loss exponents a. We find that the average ergodic 


^There are 6 cases of the allocations of spectra in three-tier heteroge¬ 
neous networks in total. However, we only analyze 4 typical cases therein. 


MU rate increases slowly with the increasing BS density 
when we fix the path loss exponent. Also, our proposed 
cell association scheme outperforms the CRE association 
scheme analyzed in na, which indicates the effectiveness 
of ours’. Path loss exponent has more effects on the average 
ergodic MU data rate when the BS density is fixed. Higher 
path loss exponent always results in higher average ergodic 
MU data rate. 

Fig. [ 3 ] compares the average ergodic MU data rate with 
respect to path loss exponent with two kinds of spectral 
allocations, i.e., (B1>B2>B3) and (B1>B3>B2). The av¬ 
erage ergodic MU data rate increases with the increasing 
path loss exponent when spectral allocation is fixed, which 
indicates that to some degree higher path loss exponent 
contributes network performance. When path loss exponent 
is fixed, (B1>B3>B2) performs better than (B1>B2>B3). 
(B1>B2>B3) represents traditional spectral allocation which 
distributes more spectral resource towards towered BSs, 
while in (B1>B3>B2) small BSs have more spectral re¬ 
source than towered BSs. It is implied that if using our 
proposed cell association scheme, small sized and high- 
density BSs should be allocated more spectral resource to 
obtain better holistic performance. 

Fig. a and fig. |5] illustrate the average ergodic MU data 
rate with respect to path loss exponent with two kinds of 
spectral allocations. We obtain similar conclusions obtained 
from fig. |3] However, the gap between the two curves 
of (B1>B2>B3) and (B2>B3>B1) in fig. 0 and the gap 
between the two curves of (B1>B2>B3) and (B3>B2>B1) 
in fig. |5] are bigger than that in fig. [3] which again indicates 
that small sized and high-density BSs should be allocated 
more spectral resource if using our proposed cell association 
scheme. 


V. Conclusions 

In this paper, motivated by the problems of existing cell 
association schemes which are merely based on one indica¬ 
tor like RSS, SINR or distance from nearby BSs, we propose 
a new cell association scheme by joint consideration of 











































Figure 4. The average ergodic MU rate with respect to path loss 
exponent with two kinds of spectral allocations in a three-tier 
heterogeneous network. 



Figure 5. The average ergodic MU rate with respect to path loss 
exponent with two kinds of spectral allocations in a three-tier 
heterogeneous network. 

SINR and traffic load in heterogeneous networks. Through 
numerical results, we find that our proposed cell association 
scheme outperforms CRE association scheme. Also, the 
results provide some insights of spectral allocation by using 
our proposed cell association scheme, which implies that 
allocating small sized and high-density BSs more spectral 
resource results in better holistic performance. 

Still, some work need to be done to further this proposed 
cell association scheme. For instance, if the number of 
candidate BSs n in each tier network is more than one, 
the corresponding analysis may be more general. And also, 
adding shadowing may make the scenario more realistic. 
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